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Abstract 

We investigate the signature of a heavy charged Higgs boson of the Minimal Super- 
symmetric Standard Model in the lepton plus multi-jet channel at the Large Hadron 
Collider with four 6-tags. The signal is the gluon-gluon fusion process gg — > tbH~, 
followed by the H~ -^ ib decay, while the main background is from gg — > tibb. We find 
that the two can be separated effectively by kinematic cuts and mass reconstruction, 
but the signal size is not very large in the end. Nonetheless, with a good 6-tagging 
efficiency, e^, ~ 50%, this channel can provide a viable signature over a limited but 
interesting range of the parameter space. 



1 Introduction 

The detection at the Large Hadron Colhder (LHC) of charged Higgs bosons would represent 
an unequivocal signal of physics beyond the Standard Model (SM). While the SM predicts 
only a neutral Higgs boson 0, any two Higgs doublet extension of it predicts a pair of 
charged Higgs bosons H^ along with three neutral ones: the CP— even H and h and the 
CP— odd A [|l|. This is true in particular for the Minimal Supersymmetric extension of the 
Standard Model (MSSM). While the SM Higgs boson may be hard to distinguish from one 
of the neutral Higgs bosons of the MSSM, the charged Higgs boson carries the unambiguous 
hallmark of the Supersymmetric (SUSY) Higgs sector. Moreover, in contrast to the case 
of super-particles, whose signature depends sensitively upon assumptions regarding the 
P— parity status and the nature of the SUSY breaking, the signature of H^ bosons is 
fairly model independent. Therefore, there has been considerable interest in looking for 
SUSY signals via the associated Higgs sector and in particular the charged Higgs bosons. 
Furthermore, the masses and couplings of all the MSSM Higgs bosons are given at tree level 
in terms of only two parameters: i.e., one of the Higgs masses (e.g., M^) and the ratio of 
the vacuum expectation values of the two Higgs doublets (tan/3). Thus the experimental 
determination of the H^ mass Mh± would go a long way to quantitatively determining the 
MSSM Higgs sector. 

In the MSSM one has a lower mass limit at tree level, Mh± > Mw±, which is not 
significantly modified by radiative corrections. There is also a comparable experimental 
limit from direct LEP2 searches for the charged Higgs boson |^. In fact, using the MSSM 
mass relations one can get a strong indirect bound on Mh± in the low tan/3 region from 
the LEP2 limit on M^, the mass of the lightest Higgs boson of the MSSM 0. As shown in 
0, however, this constraint can be evaded in modest extensions of the MSSM involving an 
additional singlet Higgs boson. Therefore it does not preclude direct searches for charged 
Higgs bosons in the low tan/3 region. 

The Tevatron and (especially) LHC hadron colliders offer the possibility to carry on the 
charged Higgs boson search to higher masses because of their higher energy reach. Here, the 
most prominent source of H^ production for Mh± < rrit (light charged Higgs boson) is top 
quark decay, t — > bH~^. The corresponding branching ratio (BR) can easily be estimated 
from the relevant part of the MSSM Lagrangian 

H^ {nih tan P ibji + nit cot P ibi), (1) 



\/2Mw± sin 9w 

written in the diagonal CKM matrix approximation. It suggests a large ibH^ Yukawa 
coupling and hence large branching fraction for the t -^ bH^ decay in the regions of low as 
well as very high tan/3, i.e., 

tan/3^1 and tan P^rrit/ mi,. (2) 



Interestingly, these two regions are favoured by 6 — r unification for a related reason: i.e., 
assuming rrih = mj. at the GUT scale one needs a large thH^ Yukawa coupling contribution 
to the Renormalisation Group Equations (RGEs) to control the rise of 771^ at low energy 
scales 1^. The dominant decay channels of light charged Higgs bosons are H^ — > cs and 
H^ — > t*h — * hbW at tan/5 ^ 1, while the H^ — > r+z/,- decay dominates for tan/3 > 1 
Pj. It may be noted here that there are already some modest limits on Mh± from the 
Tevatron top quark data in the two regions (§). The search can be extended over a wider 
region of Mh± and tan/3 at the forthcoming Tevatron upgrade (TeV-2), by exploiting the 
distinctive r polarization in H^ decay p|, ^]. Moreover the detection range can be enlarged 
to encompass practically the entire Mh± < fnt region at the LHC [^ |10|- 

We shall investigate here the prospect of charged Higgs searches at the LHC in the 
opposite case, when Mh± > rrit (heavy charged Higgs boson). The dominant decay mode 
of such a particle is H^ -^ tb, which suffers from a large QCD background. Therefore, it 
is not surprising that detecting a charged Higgs boson heavier than top has been generally 
regarded as very hardQ. The largest signal cross sections at the LHC are expected to come 
from the associated production of H"^ with top (anti) quarks, followed by its dominant decay 
mode, i.e., 

gb -^ tH' -^ tib, (3) 

and 

gg -^ tlH' -^ tibb. (4) 

After the decay of the top pair one expects 4 6— jets in (^) and 3 in (^), where the accom- 
panying b sea-quark does not take part in the hard scattering and hence escapes detection 
by being produced close to the beam direction. The charged Higgs signal from @ and (|) 
were investigated in |jT^ and [jl5|, respectively, assuming triple 6-tagging. The dominant 



(QCD) background processes in either case are 

99 -^ tibb, (5) 

99 -^ tigg + tiqq, (6) 

with one or more of the light parton jets in the latter process misidentified as b. Both 
the analyses found fairly viable H^ signals in the two above mentioned regions of tan/3, 
corresponding to a large H~^ib Yukawa coupling. Recently the associated production of 
charged Higgs and W^ bosons has been investigated in |T^. Being a second order elec- 
troweak process, however, the size of the resulting signal is smaller than those of reactions 

^Recent new insights into the problem can be found in Ref. |L1| , where the exploitation of top quark 
polarisation effects was advocated. For an alternative approach, using the much suppressed but somewhat 
cleaner H^ -^ tv decay channel, see Ref. |lj]. In the same spirit, the mode H^ -^ W^h is currently being 
investigated in Ref. flJ. 



(l^)-(^. Furthermore, it suffers from an overwhelming irreducible background induced by 

top-antitop production and decay [l^ . Therefore it does not seem to offer a useful if ^ signal 

at the LHC Similarly, the production of charged Higgs scalars in association with b quarks 

18| is burdened by a large combinatorial background as well as a formidable multi-jet QCD 



noise. 

We shall present here an analysis of H^ signals at the LHC produced via the gluon- 
gluon fusion process (|D0, assuming all four 6— quark jets to be tagged. Of course the size of 
this signal will be smaller than in the 36-tagged channel. However, we have verified, while 
computing the ^f^f-initiated process, that, by imposing a transverse momentum cut of 20 to 
30 GeV on the spectator 6— jet in reaction (^, the typical loss of signal is not dramatic: 
e.g., about a factor of 2 to 3, for Mh± ~ 300 GeV. Besides, it is clear from the analyses in 
Refs. [0, |T3| that the background (^) with double 6-mis-tagging will be very small. Indeed, 
we have explicitly confirmed this by recomputing the processes (|^)-(|^) from scratch. The 
same conclusion applies to the the case of single 6-mis-tagging in events of the type: 

gb -. tigb. (7) 

Therefore we only need to worry about the tibb background (J^). Moreover, we shall see 
below that the kinematics of the pair of 6— jets accompanying the ti pair is expected to be 
rather different for the background, as compared to the signal (^). Consequently such a 
background can be suppressed effectively by using suitable kinematic cuts on these 6— jets. 
As a result, we find a cleaner (but smaller) charged Higgs boson signal in the 46-tagged 



channel than in the 36-ones considered in [[14| , |T^. In the following Section we present 
the main steps in the calculation of the signal and background cross sections. The event 
selection strategy and choice of kinematic cuts are outlined in Section |^. The quantitative 
discussion of our analysis is presented in Section ^. Finally, we summarise the main results 
and present our conclusions in Section ^ 

2 Calculation of signal and background cross sections 

The relatively large number of Feynman diagrams involved in the signal process, gg -^ tbH^ 
(^), renders the algebraic expressions for the corresponding 'squared' matrix element rather 
long and cumbersome. This drawback of the trace method can be obviated by resorting 
to helicity amplitude techniques, which allow one to write down the 'complex' amplitudes 
in more compact form. In expressing the helicity amplitudes we have made use of the 
formalism described, e.g., in Ref. [^, to which we refer the reader for technical details. 



Here, we briefly outline the procedure and give all the formulae needed to implement our 
calculations in a numerical program, as they have not appeared in the literature previously. 



^For some early numerical studies of the on-shell 2^3 production, see Ref. pi 
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y{pi,^i;P2,\2;cR,CL) 


++ 
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+- 


ClS{ + ,Pi,P2) 



Table 1: The Y functions for the two independent hehcity combinations in terms of the 
functions S, t] and /i defined in the text. The remaining Y functions can be obtained by 
flipping the sign of the helicities and exchanging + with — in the S functions and R with 
L in the chiral coefficients. 



First, one needs to introduce some spinor functions, Y and Z [^, ^, ^], given in 
Tabs. |I|-0, which can be defined in terms of the following quantities (e°^^^ = 1 is the Levi- 
Civita tensor): 

o (Pi ■ ko){p2 ■ ki) - (pi ■ ki){p2 ■ ko) + ie^.up^k'^k'iplpl 



S{ + ,Pi,p2) 
S{-,PUP2) 



^1^2 



S{ + ,P2,Pl 



-■- ) 
Vi 



Vi 



2fe ■ ko) 



(9) 



where Pi and rrii represent the four-momentum and mass of the particle i (for which pf = mf] 
and ko and ki are two arbitrary four-vectors such that 



ko ■ ko 



0, 



ki ■ ki 



ko ■ ki = 0. 



(10) 



In the first of the two expressions in eq. (^ the signs + and — refer to particles and to 
antiparticles, respectively. 

In Tabs. |l] and ^ the (chiral) coefficients cr and c^ are those entering the fundamental 
fermion-fermion-boson vertices, described through the expressions 



and 

with 



1 + 75 „ 1-75 



(11) 

(12) 
(13) 



" 2 ' ^ 2 ' 

the chiral projectors. 

If we make the following assignments for the momenta pi (with i = 1, ...5), which we 
assume incoming in the initial state and outgoing in the final state, and helicities Aj (with 
i = 1, ...4) of the external particles in the 2 — >■ 3 reaction 



gipi, Ai) + g{p2, A2) -^ t{p3, A3) + 6(p4, A4) + H {p5 



(14) 



A1A2A3A4 


Z{pi, \i;p2, A2; P3, A3; P4, A4; cr, cl; c'j^, c^) 


+ + ++ 


-2[S{ + ,P3,Pl)S{-,P4,P2)c'jiCR - l^ll22r]3ri4c'^CL - r/ir/2/i3/i4CiCi?] 


+ + +- 


-2ri2CR[S{ + ,p4,Pl)l^3c'L - S{ + ,P3,Pl)fi4c'ji] 


+ + -+ 


-2r]iCR[S{-,P2,P3)l^4c'L - S{-,P2,P4)I^3C'r] 


+ -++ 


-2ri4C'ji[S{ + ,P3,Pl)l^2CR - S{ + ,P3,P2)I^iCl] 


+ + — 


-2[S{ + ,Pi,P4)S{-,P2,P3)c'lCr - ^il^i2V3V4C'LCL " ^1^2/^3/^4C^Ci^] 


+ -+- 





+ --+ 


-2[^l^4r]2r]3c'LCL + /i2/i3'7l^4 C^Ci? - /i2/i4^1?73CiCi? - IJ,iIJ,3V2V4c'rCl] 


+ 


-2r/3c'i[S( + ,p2,P4)/ilCL - S{ + ,Pi,P4)H2Cr] 



Table 2: The Z functions for all independent helicity combinations in terms of the functions 
S, 7] and /x defined in the text. The remaining Z functions can be obtained by flipping the 
sign of the helicities and exchanging + with — in the S functions and R with L in the chiral 
coefficients. 



so that p\ = pI = 0, pI = mf (with Ts = Ft ^ 0), pi = ml (with r4 = F;, = 0) and 
pI = M]j± , then the Feynman amplitudes T/ can be written (apart from a phase, a factor 
g'^e and neglecting the colour matrices) as (here and below, {A} refers cumulatively to the 
hehcities of the external particles Aj with i = 1, ...4, and J2{x} to the summation over all 
their possible combinations) 



n 



{A} 




(15) 



+ N^N2D4{P3 + P5)D4{P2-P4) E E^AEE 

4=1,2,4 j=2,4 A=±A'=± 

^(P3, A3; Pi, A; cf,cf) 
X Z{pi, A; pj, A'; pi, Ai; gi, Ai; 4, 4; 1, 1) 
X Z{pj, A'; P4, -A4; P2, A2; g2, A2; cfj, ci; 1, 1), 



T. 



{A} 




(16) 



- N,N2D,{p,-p,)D,{p2-p,)Y. T.brbjY. E 

1=1,3 j=2,4 A=±A'=± 

Zips, A3; Pi, A; pi, Ai; qi, Ai; cfj, c^; 1, 1) 
X F(p„A;pj,A';cf,cf) 
X Z{pj, A'; p4, -A4; P2, A2; g2, A2; cfj, c[; 1, 1), 



T. 



{A} 




(17) 



+ N,N2D,ips-pi)Dsip,+P5)T. E M, EE 

1=1,3 i=l,2,3 A=±A'=± 

ZiP3, h;Pi, A; Pi, Ai; qi, Ai; 4, c|; 1, 1) 
X Z{pi, A; pj. A'; pa, A2; ^2, A2; 4, 4; 1, 1) 
X Y{pj,\';pi,-\i;c^,cf), 



z 
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t/^^(i 



2), 



T. 



{A} 
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t}'\i ^ 2), 



(19) 



T, 



{A} 




.{A} 



Tr(1^2), 



(20) 



onro^ 



TJ^^ = ^%Tnn5X """ = (21) 



j=l,2,4 A=±A'=± 

>"(P3, A3; Pi, A; cf,cf) 
X [Y{p„X■p^■X'■l,l)Y{p^,X'■p,,-X,■c'^,ci) 
xZ{pi, Ai; gi, Ai; pa, A2; ^2, A2; 1, 1; 1, 1) 
+2 Z{p2, A2, g2, A2; Pi, A; P4, -A4; 1, 1; c% c[) 
xr(pi,Ai;p2;A';l,l)r(p2,A';gi,Ai;l,l) 
— same(l ^^ 2)], 



Tf^ = >^nnr/__^ = (22) 

2 "4 

- mN2D,{p^+p,)D{p,+p2) T. biY.Il 

1=1,2,3 A=±A'=± 

[Y{P3, A3; Pi; A'; 1, l)Y{pi, A'; Pi, A; cf^, c[) 
xZ{pi, Ai; gi, Ai; pa, A2; ^2, A2; 1, 1; 1, 1) 
+2 Z(p2, A2, g2, A2; P3, A3; Pi, A; 1, 1; cfj, 4) 
xF(pi,Ai;p2;A';l,l)y(p2,A';gi,Ai;l,l) 
— same(l <-> 2)] 
X >"fe,A;p4, -A4;c^,cf), 

where we have introduced the coefficients 61 = 62 = — ^s = — &4 = 1, to distinguish incoming 
and outgoing particles, and the propagators Di{p) = l/ijP' — mf + irrii'yi) (with 7^ = Fj 
if p^ > and 7i = otherwise^) and D{p) = 1/p^. In eqs. (P^)-(P^), qi (with i = 1,2) 
is an arbitrary four- momentum not parallel to Pi (i.e., q^ 7^ apj with a constant) and 
Ni = [4(gj -pi)]""^^^; see Ref. [^ for more details. 

•^That is, we include a finite quark width only in resonant propagators. 



The coefficients cr and cl for the vertices relevant to such a process are 

\ 

(4'Cf) = — 7=— r^("^fetan/3,mtCot/5) (23) 

for the charged Higgs, see eq. ([l|), and simply 

(4,4) = (1,1) (24) 

for the gluon to quarks couplings. As usual, 6w represents the Weinberg angle. Here, 
the bottom and top quark masses entering the Yukawa interaction are those defined at 
the propagator pole, i.e., rrib^t = T^b,t{jnb,t) , where fhb^t{Q) are the running masses at the 
(energy) scale Q (see below). Also note that in eqs. (|23|)-(^4D we have factored out the 
overall couplings — ie and —igs of the Lagrangian. 

As for the colour structure of our process, one should notice that in this case there are two 
basic combinations of colour matrices associated with the Feynman graphs ([r5|)-(p2D, that 
is (tH^)i^i^ and (t^t^)i^i^, where A^i^) and -B(i4) identify the colours of the gluons(quarks) 
1(3) and 2(4), respectively. In fact, it should be recalled that the colour terms associated 
with the triple-gluon vertices are nothing but the structure constants f^^^ of the SU{Nc) 
gauge group, for which 

[t yt \i3i4 = [t t ji3J4 ~ [t t )i^i^ = 1/ ''J3i4, I20j 

X being in our case the colour label of the virtual gluon. Therefore, one can conveniently 
group the original eight Feynman amplitudes as follows 

Mf^ = J2t!^^ + J2tI^^ (26) 

1=1 i=7 

1=4 i=7 

The standard form of the matrix element (ME) squared, summed/averaged over the fi- 
nal/initial spins and colours, is then 

IM^gg - fbH-) = ^ E E Mf >Mf >*a„ (27) 

where Cij is a 2 x 2 colour matrix with elements 






C++ = C'__ = ^(Ar^-2 + ^) (2^ 



c = c -^^r ^ 



Nc = 3 being the number of colours in QCD. Finally, notice that the ME for the charge 
conjugated process 

9{Pu Ai) + 9{P2, A2) ^ KP3, A3) + t(p4, A4) + H+{p5), (29) 

now with pI = ml and pi = m^ (and, correspondingly, Fs = and F4 = F^), can be obtained 
by the simple replacement c^ ■«->^ c|f . 

The above formulae refer to the 2^3 process of on-shell H^— and t— production. 
In reality, both these particles eventually decay inside the detectors, so that one ought 
also to consider their decay signatures. We have included these decays by convoluting 
the 2^3 (unpolarised) production ME with the 3— and 4— body decay MEs of top and 
Higgs, respectively, which are well known and can be found in the literature. In doing 
so, we introduce several simplifications^. Firstly, we neglect spin effects in the top decays. 
Secondly, we do not consider interference effects between diagrams involving H^ and H~ 
production. Thirdly, Fermi-Dirac interferences due to indistinguishibility of 6— quarks (or, 
equivalently, of 6— antiquarks) in the final state are not included. However, while noticeably 
simplifying the numerical evaluation, we have checked that these approximations do not 
spoil the validity of our analysis. In fact, we have also produced the exact 2 — > 8 ME for the 
signal process, including all the above spin and interference effects, by means of the HELAS 
[p4| subroutines and the MadGraph ||2^ package, and compared its yield to that of the 
simplified implementation. We have always found a good agreement between the two, with 
residual effects surviving only in differential spectra which are of no concern in our analysis. 
Indeed, for the typical quantities we shall investigate (transverse momenta, pseudorapidity, 
multi-jet invariant masses, etc.), the results generally coincide within numerical accuracy. 

The HELAS libraries and MadGraph have also been used to generate the background 
process (5a), as any analytical expression for the latter would be much too cumbersome. 
Indeed, thirty-six different Feynman diagrams are involved (actually twice that if one con- 
siders also those induced by the above mentioned Fermi-Dirac statistics), with ten external 
particles. 



Both signal and background MEs have been integrated by means of VEGAS [^, with a 
careful mapping of the phase space, to account for the various resonances. In some cases, 
where the multi-dimensional integrations (21 in total, for the 2 — >■ 8 implementation) are 
more problematic, the VEGAS results have been cross-checked against those obtained by 
using RAMBO [^ as well as the DOIGCF and DOIGZF NAGLIB subroutines. 

In addition to the phase space integration, one also has to convolute in the 
{x, (5^)-dependent Parton Distribution Functions (PDFs) for the the two incoming glu- 
ons. These have been evaluated at leading-order, by means of the packages MRS- 
LO(05A,09A,10A,01A,07A) ^. The Q^ used for the latter was the CM energy (squared) 

^Note that in the numerical simulations we allow the top quark and Higgs bosons to be 'off-shell', that 
is, Pgr^i ^ m^ and p^ ^ Mfj± in the above formulae for process (p^[(p9|)]. 



at the parton level, i.e., s = X1X2S, with ^/s = 14 TeV taken as the CM energy for the 
LHC The same choice has been adopted for the scale of the strong coupling constant a^, 
evaluated again at lowest order, with Aq^^ chosen in accordance with that of the PDF set 
used. 

In addition, notice that before the signal and background cross sections can be computed 
reliably, one must take into account the effects of higher-order QCD corrections to the tree- 
level processes (see the discussion in Ref. |]15|). The full next-to-leading order corrections are 
as yet unknown for the processes we consider. For the case of the Higgs signal, the dominant 
effects can however by mimicked by adopting the pole masses in the Higgs-fermion coupling 
entering the production process, see eq. (]23|), rather than the running ones ^^. As for the 
tibb background, we estimate their effect by including an overall i^-factor of 1.5 throughout 
this study []T3|. 



= rris 


= mc = 0, 


nit = 


175 GeV, 


Tz = 


= 2.5 GeV, 


Ty/ = 


= 2.08 GeV 



The numerical values of the SM parameters are (pole masses are assumed): 

mi = rriy^ = rriu = rrid 

mb = 4.25 GeV, 

Mz = 91.19 GeV, 

Mw± = 80.23 GeV, Tw = 2.08 GeV. (30) 

For the top width F^, we have used the LO value calculated within the MSSM (i.e., F^ = 1.55 
GeV for Mh± ^ rrit). Furthermore, we have adopted Ma and tan/5 as independent input 
parameters defining the Higgs sector of the MSSM at LO. The charged Higgs width Th± has 
been computed by means of the program HDECAY, which indeed requires M4, rather than 
Mh± , as mass input [Q (the masses of the super-partners of the ordinary matter particles 
were fixed well above 1 TeV, so they enter neither the top nor the charged Higgs decay 
chain as real objects and render the virtual SUSY corrections to the H^tb vertex negligible 
pT|). As this program uses running quark masses (in the modified Minimal Subtraction 
scheme) in evaluating the decay width of the H^ — > tb channel we have, for consistency, 
used running values for all such quantities in the H^tb couplings entering our decay MEs 
for the signal (but not in the propagators and in the phase space, for which the pole masses 
nib = tfibijni,) and rrit = fhtirrit) of eq. (|30|) have been used). 

Finally, notice that we stop our analyses at the parton level, without considering frag- 
mentation and hadronisation effects. Thus, jets are identified with the partons from which 
they originate and all cuts are applied directly to the latter. In particular, when selecting 
6— jets, a vertex tagging is implied, with finite efficiency e?, per tag. As four 6— jets will be 
required to be tagged, the overall efficiency will be e^, by which both signal and background 
rates will eventually have to be multiplied. For simplicity, we assume no correlations be- 
tween the four tags, nor do we include mis-identification of light-quark (including c— quark) 
jets produced in the W^ decay as 6— jets. 

10 



3 Selection strategy 

In this Section we describe the kinematic procedure adopted to disentangle charged Higgs 
events (^ from the background (^ in the tibb channel. First, one of the top quarks is 
required to decay leptonically (t — > blu) to provide a hard lepton (£ = e, /i) trigger and avoid 
the QCD background, while the other decays hadronically {i -^ bqq'), with {q,q' 7^ b,t). 
The resulting branching fraction is 2 x 2/9 x 2/3, with the factor of two accounting for the 
fact that each of the PF^'s can decay either leptonically or hadronically. We assume that 
the charge of the 6— jet is not measured. Thus, the signature we are discussing is effectively 

46 + 2 jets + £± + pLs, (31) 

where the 2 un-tagged jets and the ^+Pmiss(= Pu) arise from the W^W~ boson pair produced 
in the ti decay. 

We note that all the decay products of the top quarks in the signal are expected to be 
hard while one of the accompanying 6— quarks (or both depending on the Mh± — rrit mass 
difference) could be soft. Therefore we impose a relatively demanding transverse momentum 
cut on the two softer 6— jets, 

pI^,^ > 20 GeV. (32) 

For simplicity, we do the same for the rest, i.e., 

pI,u,jMM > 20 GeV. (33) 

However, we will present some results also for the case of a 30 GeV cut in transverse 
momentum (on both jets and leptons), since the latter threshold is believed optimal in 
increasing the 6-tagging efficiency at high luminosity (this is needed in order to render the 



46-signal of a charged Higgs boson statistically significant): see Ref. [^. Furthermore, we 
require the pseudorapidity of jets and leptons to be 

\7]hj^i±\ < 2.5, (34) 

and allow for their detection as separate objects by imposing the following isolation criteria: 

^Rbb,bj,jj,bl^,ji± > 0.4, (35) 

by means of the variable 



Ai?,, = ^(Ar/,,•)2 + (A0,,.)^ (36) 

defined in terms of relative differences in pseudorapidity rjij and azimuth (pij, with i j^ j = 

We simulate calorimeter resolution by a Gaussian smearing of p^ (without shower 
spreading and with uniform pseudorapidity/azimuth segmentation), with (o"(p^)/p"^)^ = 

11 



(0.6/yy^)2 + (0.04)2 fo^ g^ii ^j^g jg^g g^^^ ^^^pT)/pT)2 _ (0.12/yy^)2 + (0.01)2 for ^i^g igp_ 
tons |]n|. The corresponding p^iss is evaluated from the vector sum of the jet and lepton 
transverse momenta after resolution smearing. 

To improve the signal/background ratio and to estimate the if ^ mass we follow a strategy 
similar to the one in Ref. |]H[ . 

(a) The invariant mass of the two un-tagged jets is required to be consistent with M^yi, 

\Mjj-Mw±\<l5GeV. (37) 

(b) The neutrino momentum is reconstructed by equating p^ = Pj^iss ^^^ fixing the lon- 
gitudinal component p^ via the invariant mass constraint M{(!,v) = My/±. The latter 
gives two solutions. If they are complex we discard the imaginary parts and they 
coalesce; otherwise both the solutions are retained. 

(c) The invariant mass formed by combining the un-tagged jet pair with one of the four 
6— jets is required to match rrit, 

|Mjjb-"^t| < 25 GeV. (38) 

If several 6— jets satisfy this constraint, the one giving the best agreement with mt is 
selected. 

(d) The invariant mass formed by combining £ and v with one of the 3 remaining 6— jets 
is also required to match rrit within ±25 GeV. If several 6— jets satisfy this, the one 
giving the best agreement with mt is selected along with the corresponding value of 

(e) The remaining pair of 6— jets may be looked upon as the hb pair accompanying the 
ti in the signal (|[) and background (|^). Note that one of these 6— jets is expected 
to come from the H"^ decay in the signal, while for the background they both come 
from a gluon splitting. Consequently, in the latter case one supposes the bb pair to 
have a smaller invariant mass. Furthermore, one may also expect the energy and the 
relative angle of the two 6's to be rather different, between signal and background. 
We compare the signal and background cross sections against M;,;,, Eb^, Eb^ (with the 
labels 1[2] referring to the more [less] energetic of the two 6— quarks) and cos6'ftfe and 
suppress the latter by suitable cuts in one or more such quantities. 

(f) Finally we combine each of the reconstructed t— quarks with each of the remaining 
6— jets to obtain 4 entries for the bt invariant mass Mf,t- For each signal point, one 
of these entries will correspond to the parent H^ mass while the other three will 
represent the combinatorial background. We plot the signal and background cross 
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sections against this quantity. The former will show the resonant peak at Mii± sitting 
on top of a broad combinatorial background, while the latter will show only a broad 
distribution in Mfe^. As we shall see below, the Breit-Wigner peak itself can help to 
improve the signal/background ratio further as well as to determine the H^ mass. 

(g) For Mh± ^ Tnt, one of the above mentioned 6— jets (i.e., the one coming from the 
H^ -^ th decay) would generally be much harder than the other. In this case, one can 
expect to reduce the combinatorial background by combining each of the top quark 
pair with the harder of the two accompanying 6— jets. This would give two values of 
invariant mass M^j for each signal point, one of which corresponds to the parent if ^ 
mass. Therefore we shall also show the signal and background cross sections against 
this quantity for Mh± > 300 GeV. 

4 Results and discussion 

Both signal (^ and background (^ cross sections are finite over the entire phase space, 
provided the 6— quark mass is retained in the calculation. Thus, as a preliminary exercise, 
we look at the total production rates of the above processes with no cuts whatsoever, as 
they would appear in the decay channel (|3l| ) to an ideal detector. This is done in the top 



plot of Fig. |l|. Here, the rates have been obtained by multiplying the 2^3 cross section 
times the relevant BRs of top and Higgs decays, thus neglecting finite width effects. The 
signal rates depend on both Mh± and tan/3, and so they are plotted as a function of the 
former for two values of the latter in the favourable regions (^. In contrast, the background 
rates are independent of both and are indicated by the arrow next to the y— axis. The insert 
in the top plot of Fig. |I| enlarges the region around the threshold region Mh± ~ fnt^ where 
the BR(i/+ -^ th) increases rapidly. A striking feature of the production cross sections 
is the apparently poor signal-to-background ratio, over all the Mh± spectrum considered, 
irrespective of tan (3. Note however that no dedicated treatment of the final state kinematics 
has yet been performed. Indeed, the number of heavy charged Higgs events produced is 
sizable up to around 800 GeV, where the total cross section at both tan/3 values is still 
around several femtobarns. At its maximum, just after the opening of the H^ — > th decay 
threshold, it can be larger by about two orders of magnitude. 

Fig. |T|has been produced by using the MRS-LO(05A) PDF set. In Tab. ^we study the 
dependence of both signal and background rates on the choice of the structure functions, 
using the other four fits of the 1998 Martin-Roberts-Stirling- Thorne LO package]^. For ref- 
erence, we have used the value tan/5 = 30, whereas five different choices of Higgs masses 
in the heavy range have been adopted. Differences in the signal cross sections are found 

^The additional fits correspond to varying the large-x gluon and a^ about their central preferred values. 
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to be within the ±25% range, with a somewhat smaller range for the background. Fur- 
thermore, changing the renormalisation (Q) and factorisation (/i) scales from their common 
default value Vs to, e.g., the sum of the rest masses in the 2^3 and 2^4 production 
process @) and (H), induces variations in the results of the same order as above. As a 
consequence, an overall error of, say, approximately 30-35% should be taken as an estimate 
of the uncertainties related to the PDFs and a^ throughout the paper. 

As the next step of our analysis, we implement the acceptance cuts (p2D -(p5D and the 
selection cuts described through steps (a)-(d) in the previous Section. The total signal and 
background cross sections after such constraints have been enforced can be found in the 
bottom plot of Fig. |l|. Note that the kinematic procedure outlined above has been helpful 
in increasing the signal-to-background ratio over all the Higgs mass spectrum (compare 
to the curves above). The signal rates have been depleted too, mainly by the p^ cut 
(|5^), dropping to a few femtobarns for values of tan/3 at the upper and lower end of the 
parameter range and Higgs masses below 700 GeV or so. For intermediate values of tan/?, 
e.g., around the minimum of the production cross section occurring at tan/3 ^ 7, prospects 
are more gloomy. In fact, the signal rates are always below 1 fb or so in this case, for any 
Mh± value. Furthermore, as here (and in the following as well) we have retained a finite 
value for Th, the reader may appreciate - by comparing this plot to the one above - the 
suppressing (enhancing) effects of a larger (smaller) Higgs width at low(high) Mh± values for 
tan/? = 40, with respect to the case tan/3 = 1.5 (the effects of a F^ 7^ are the same in 
all cases). As intimated in the previous Section, we also have considered the case of a 30 
GeV cut in all transverse momenta (including the missing one). The insert in the bottom 
plot of Fig. |l] presents the consequent suppression on the signal rates for, e.g., tan/? = 40. 
Far from the Mh± ^ mt + nib threshold, the ratio between the two cross sections stabilises 
at just below three (irrespectively of tan/3). For the background, the reduction is slightly 
higher, a factor of four or so. 

Continuing with step (e) of our plan, we next investigate the mass, angular and energy 
behaviours of the two 6— quarks accompanying the ti pair, after the acceptance and selection 
cuts have been implemented. The relevant plots can be found in Fig. |^. For reference, 
the value chosen for tan/? is 40, whereas for the charged Higgs masses we have taken 
Mh± = 214(310) [407]{506} GeV, corresponding to M^ = 200(300) [400]{500} GeV. (Spectra 
look rather similar if a pj± ^j, 6 > 30 GeV cut is enforced instead.) As already foreseen, one 
can appreciate significant differences in all four quantities considered, as long as Mfj± is well 
above m,j. If Mh± is not much larger than 200 GeV, the signal and background distributions 
are similar except for the angular variable. In contrast, if Mh± > 300 GeV, one can see a 
greater discriminatory power in each of these variables. 

In order to enhance the signal-to-background ratio, especially in the very heavy Higgs 
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mass region, we therefore adopt the following additional constraints on the two 6— jet system: 
Mbb > 120 GeV, cos9bb < 0.75, Ef,, > 120 GeV. (39) 

The resulting signal cross sections are shown in Fig. ^ for tan/3 = 1.5 and 40. along with 
those of the background. The signal-to-background ratio has increased, but the background 
remains larger. However, notice the very little loss of signal at very large Higgs masses. In 
fact, the Higgs rates remain above 1 fb for Mh± ~ 800 GeV. For a transverse momentum cut 
of 30 GeV throughout, the typical suppression on the signal (away from threshold) is again 
about a factor of 3 (see the insert in the figure), now similar to the case of the background. 

To enhance the relative rates further and estimate the charged Higgs boson mass, we 
reconstruct the M^t invariant mass by combining each of the reconstructed top quarks with 
each of the accompanying 6— jets, as described in step (f), a quantity that we label M^^H). 
The resulting spectra are presented in the top plot of Fig. H for tan /? = 40 and three selected 
values of the H^ mass along with the background. The signal distribution clearly shows a 
resonance at Mi{H) k, Mh± sitting over a combinatorial continuum, while the background 
spectrum is broader and tends to concentrate at values of M^i^H) below 300 GeV. One can 
sharpen the resonances at the cost of reducing the size of the signal by combining each of the 
reconstructed t— quarks with the harder 6— jet, as in step (g) (the corresponding invariant 
mass is labeled as M2{H)). The distributions that we obtain in this way are given in the 
bottom part of Fig. |. 

One sees from the figure that the signal-to-background ratio is (much) greater than one 
in the neighborhood of the resonant peaks. This is an advantage of the 46— tagged channel 
over the 36— tagged case considered in Ref. 11^, where the backgrounds were found to exceed 



the signal. In contrast, the size of the Higgs cross section is smaller in our case, because 
of the kinematic suppression induced by requiring the detection of the fourth 6— quark. 
However, with an annual luminosity of 100 fb~^, expected from the high luminosity option 
of the LHC, and a very good 6-tagging efficiency, one would obtain a clearly viable signal. 
To illustrate this we show the signal rates on the right-hand scale of Fig. ^, for an optimistic 
6-tagging factor of e^ = 0.1, corresponding to e?, = 56%. Such a high value is now considered 
realistic for the TeV-2 run at FNAL and can be achieved by combining the silicon vertex 



and the lepton tagging efficiencies |]33[. It is also close to the optimistic expectation of a 



6-tagging efficiency of about 50% even for the high luminosity run of the LHC |3^ . 



For the above efficiency and luminosity, one would obtain between ten and hundred Higgs 
events per year with signal-to-background ratios above one for M^j- as large as 800 GeV (for 
tan/3 = 40), as shown in Tab. ^. This shows the predicted number of events in a window of 
80 GeV centered around the Higgs resonances for both signal and background, together with 
the corresponding statistical factors S/\fB. By looking at those rates, one would expect an 
accessible signal for Mh± ^ 600 GeV in the high tan/? region ( ^ 40). Similar results also 
hold for the case of low tan/3 values ( ^ 1.5). 
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However, given the not too large rates of the surviving Higgs events, the actual size of 
the MSSM parameter space that can be covered strongly depends on the 6-tagging efficiency, 
eb. For instance, changing it from 56% to 40% would result in a reduction of ef by a factor 
of four. This corresponds to a suppression of the S/yB rates of Tab. ^ by a factor of two. 
A similar effect would occur if a 30 GeV cut in transverse momenta of both missing and 
observable particles is enforced, as opposed to the 20 GeV value advocated here. In this 
case, the suppression would be somewhat smaller though, about a factor of 3 on the event 
rates and 1.7 on the statistical significances. 

Before concluding, we would like to come back to and justify what we have mentioned 
in the Introduction: that the size of the backgrounds (§)-(0) is generally smaller than 
that of process (||) considered so far, for the 6-tagging efficiencies assumed in this paper. 
Rather than re-running all the simulations for each of these additional final states, we have 
assumed the two (anti)top quarks to have already been reconstructed, with similar efficiency 
in each case. This way, we can compute the 2 — > 4 cross sections (H)-® and apply the 
transverse momentum, pseudorapidity and separation cuts of Section ^to the jet-jet system 
accompanying the ti pair, alongside those of eq. (^). We do so also for the case of process 



(1^). (For all such computations we have resorted again to MadGraph ||25|.) By adopting 
eb = 0.56 and assuming e^, q-Lb = 0.01 to be the rejection factor against non-6— jets, then the 
rates obtained this way scale as follows: 

a{gg^tibb) : a{gb ^ tigb) : a{gg ^ tigg + tiqq) ^ 3.38 : 0.56 : 0.48. (40) 

Thus, the singly and doubly 6-mis-tagged backgrounds are both one order of magnitude 
smaller than the pure 46-process. The same applies also to the other channel contributing 
to a possible double 6-mis-tagging, i.e., 

gq -^ tigq, (41) 

where q ^ b, whose production rates are in fact comparable to those of process (|^) (about 
25% smaller). Altogether, they would add a 40% or so contribution to the tthb background. 
Some of these backgrounds are likely to be enhanced by a larger probability of a c— quark 
jet faking a 6-tag, which we have not taken into account. However, their inclusion would 
not change our main results, so that we have left these reactions aside for the time being. 

5 Summary and conclusions 

The discovery at the LHC of charged scalar particles would definitely confirm the existence 
of new physics beyond the SM. In this respect, a very special role is played by the charged 
Higgs bosons of 2-Higgs Doublet Models (2HDMs), as their production and decay dynamics 
can entirely be described at tree level by only two parameters. However, the feasibility 
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of their detection at the LHC has always been far from certain if the mass of the new 
particles is much larger than the top mass. Therefore, a high-on-the-list priority is to devise 
phenomenological strategies that would allow one to meet the difficult challenge of their 
detection at the LHC collider. 

We have contributed here to this task by considering the production and decay of charged 
Higgs scalars of the MSSM in the channel gg — > tbH~ + ibH^ -^ bbti -^ bbbbW^W~ , in 
which one W^ decays hadronically and the other leptonically. The major feature of our 
analysis, as compared to others carried out in the past, is that all four 6— quarks present 
in the final state are required to be recognised as such. The advantage of this procedure is 
that it allows one to exploit the differences existing between signal and background in the 
kinematics of the heavy quark jets. In fact, in the dominant background, gg -^ tibb, the 
two 6— quarks produced in association with the ti pair are soft, collinear and at rather low 
invariant mass. In contrast, in the Higgs signal, at least one of the two is expected to be 
energetic and isolated, as long as Mh± is significantly larger than rrit. The disadvantage is 
that the additional 6— quark produced in Higgs events has rather low transverse momentum, 
so that its detection imposes a sizeable loss of signal. 

By exploiting a selection procedure that allows one to reconstruct both top and antitop 
masses and after imposing tight constraints on the two 6— quarks accompanying the ti pair, 
we do see Higgs peaks appearing on top of a flat combinatorial background and also above 
the continuum from tibb events. Their statistical significance is such that viable signals can 
be obtained for charged Higgs boson masses up to 600 GeV or so, when tan/? is either below 
1.5 or above 40, with a total number of Higgs events of the order of a several tens every 100 
inverse femtobarns of luminosity. 

Such mass coverage is significantly higher than that achieved in previous analyses based 
on a triple 6-tagging. However, we should stress that these conclusions rely on a high, 
though not unrealistic, 6-tagging efficiency, of about 50% per single 6— jet, and a transverse 
momentum cut on jets and leptons, of 20 GeV, somewhat lower than the threshold normally 
considered. If such performances can be achieved by the LHC detectors while the machine 
is running at high luminosity, then the 'lepton plus 46' channel represents a profitable new 
means to access such elusive yet crucial particles over significant portions of the MSSM pa- 
rameter space. Besides, this channel will be very useful for the measurement of the charged 
Higgs boson parameters, given the higher purity of the signal in this case. Certainly, our re- 
sults are sufficiently optimistic to justify a more detailed detector-level study, incorporating 
hadronisation of the final state partons, the effects of jet identification and reconstruction 
as well as a full background simulation. 
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Signal, tan/5 = 30 (fb) 


Ma{h±) 


05A 


09A 


lOA 


OlA 


07A 


200(214) 
300(310) 
400(407) 
500(506) 
600(605) 


217.14 ±0.18 
123.20 ±0.10 
60.414 ±0.059 
30.802 ± 0.044 
16.341 ±0.029 


226.67 ±0.20 
130.79 ±0.11 
65.256 ±0.066 
33.895 ±0.047 
18.268 ±0.032 


199.11 ±0.15 
110.90 ±0.086 
53.410 ± 0.049 
26.742 ± 0.040 
13.977 ±0.026 


204.48 ±0.17 
116.42 ±0.10 
57.333 ± 0.057 
29.481 ± 0.043 
15.670 ±0.029 


222.79 ±0.18 
125.76 ±0.11 
61.360 ±0.059 
31.252 ±0.045 
16.455 ±0.029 


Background (fb) 




05A 


09A 


lOA 


OlA 


07A 




1863 ± 56 


1949 ± 59 


1706 ± 51 


1755 ± 53 


1901 ±57 


MRS-LO [Q = ll=^/I] 


Ab + 2 jets + i^ + p^i33 No cuts 



Table 3: Production cross section for the signal (^) (and its charge conjugate) in the decay 
channel ( pi]) (top), for tan/3 = 30 and five values of M^(j|^±) (given in GeV) in the heavy 
mass range, as obtained by using five different sets of PDFs. Corresponding rates for 
the background (|5|) yielding the same signature (^Tj) are also given. No cuts have been 
implemented. The renormalisation and factorisation scales are set equal to the partonic 
CM energy. Errors are as given by VEGAS. 



Number of events per year 


Mh± ± 40 GeV 


S 


B 


S/^ 


310 


127.80 


105.14 




12.46 




57.80 


41.90 




8.92 


407 


88.56 


67.63 




10.76 




53.78 


39.21 




8.58 


506 


51.46 


38.74 




8.26 




36.32 


26.86 




7.00 


605 


29.43 


21.49 




6.34 




22.70 


16.58 




5.57 


704 


17.09 


11.98 




4.93 




14.02 


9.89 




4.45 


803 


10.03 


6.75 




3.86 




8.61 


5.81 




3.57 


MRS-LO [Q = ii = ^s] 


46 + 2 jets + £^ + p^iss After all cuts 



Table 4: Number of events from signal @ (and its charge conjugate) {S) and background 
d^) {B) in the decay channel (^Tj), along with the statistical significance {S/\fB), per 100 
inverse femtobarns of integrated luminosity, in a window of 80 GeV around a few selected 
values of Mu± (given in GeV) in the heavy mass range, with tan/3 = 40. Four 6— jets are 
assumed to be tagged with overall efficiency el = 0.1, i.e., eb = 56%. All cuts discussed 
in the text, (P^-(P3D, (a)-(d) and (^), have been implemented. The renormalisation and 
factorisation scales are set equal to the partonic CM energy. The first row corresponds to 
the Mi{H) distribution whereas the second refers to the M2{H) one (see Fig. H). 



b t H + c.c. ^bbtt^4b2jl + rnisKing energy/nionienti: 



10* 



10 



3 _ 



10' 



2 _L 




200 



300 



500 



600 



400 
Mh± (GeV) 

b t H + c.c. ^bbtt^4b2jl + missing energy/momentuin 




MRS-LO(05A) [ Q = /i = V s ] 



200 



300 



400 
Mh± (GeV) 



500 



600 



Figure 1: Production cross section for process @ (and its charge conjugate) in the decay 
channel (|3TD as a function of Mh± in the heavy mass range, for some discrete values of tan/3, 
for the case in which no (upper figure) and acceptance plus selection (lower figure) cuts have 
been implemented. In the insert of the upper plot, we enlarge the rates in the vicinity of 
Mh± = TTLt. In the insert of the lower plot, we present the ratio between the above signal 
cross section for tan/5 = 40 and the corresponding one obtained by adopting a threshold 
of 30 GeV in (il)-(|3|). The PDF set used was MRS-LO(05A) with renormalisation and 
factorisation scales set equal to the partonic CM energy. The arrow represents the size of 
the background (|^) yielding the same signature (|3l|). No 6-tagging efficiency is included. 
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Figure 2: Differential distributions in invariant mass (top-left), in cosine of the relative 
angle (top-right) and in energy of the most (bottom-left) and least (bottom-right) energetic 
of the two h quarks accompanying the tt pair, for process (HI) (and its charge conjugate) 



in the decay channel (0), for four selected values of Mii± in the heavy mass range, with 
tan/3 = 40. Acceptance and selection cuts have been implemented here. The PDF set used 
was MRS-LO(05A) with renormalisation and factorisation scales set equal to the partonic 
CM energy. The fifth (dot-dashed) curve represents the shape of the background (^ yielding 
the same signature (|3T|) . All distributions are normalised to unity. No 6-tagging efficiency 
is included. 
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Figure 3: Production cross section for process @ (and its charge conjugate) in the decay 
channel (|3T|) as a function of Mh± in the heavy mass range, for two discrete values of 
tan/3. Acceptance and selection cuts have been implemented here, along with the additional 
cuts ( p9| ) on the 26— system accompanying the ti pair. In the insert, we present the ratio 
between the above signal cross section for tan/? = 40 and the corresponding one obtained 
by adopting a threshold of 30 GeV in (§g)-(|S|). The PDF set used was MRS-LO(05A) with 
renormalisation and factorisation scales set equal to the partonic CM energy. The arrow 
represents the size of the background (|^) yielding the same signature (^Tj). No 6-tagging 
efficiency is included. 



b t H + c.c. ^bbtt^4b2jl + missing energy/momentum 



■d 

2 



0.30 



0.25 — 



0.05 W 



0.00 




MRS-L0(05A) [ Q = fj. = -J s 

Mh+ = 310 GeV 

Mh+ = 407 GeV 

Mh+ = 506 GeV 

Mh+ = 605 GeV 

Mh+ = 704 GoV 

Mh+ = 803 GeV 

tan/3 = 40. 



— 0.020 



200 



400 



800 



1000 



0.000 



600 
M4(H) (GeV) 

b t H + c.c. ^bbtt^4b2jl + missing energy/momentum 

0.150 I , , , , , , , , , , , , , , , , , , , , , m 0,015 



0.125 




200 



400 



MRS-LO(05A) [ Q = /i = V s ] 

M„+ = 310 GeV 

M„+ = 407 GeV 

M„+ = 506 GeV 

M„+ = 605 GeV 

M„+ = 704 GeV 

M„+ = 803 GeV 

'( tan/3 = 40. 

I 
I I 

I ' /■■■. 

600 800 1000 

V-A^) (GeV) 



0,010 



Figure 4: Differential distributions in the reconstructed charged Higgs mass for process (D 
(and its charge conjugate) in the decay channel ( pT] ) for six selected values of Mh± in the 
heavy mass range, for tan/? = 40. Acceptance and selection cuts have been implemented 
here, along with the additional cuts ( P^ on the 26— system accompanying the ti pair. The 
PDF set used was MRS-LO(05A) with renormalisation and factorisation scales set equal 
to the partonic CM energy. The seventh (dot-dashed) curve represents the shape of the 
background (|^) yielding the same signature (0). Normalisation is to the total cross sections 
times the number of possible '26 + 2 jet mass' combinations: four (top) and two (bottom). 
The right-hand scale corresponds to a 6-tagging efficiency factor of ef = 0.1, i.e., eb = 56%. 



